OPTIMAL ESTIMATES FROM BELOW FOR BIHARMONIC 
GREEN FUNCTIONS 
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Abstract. Optimal pointwise estimates are derived for the biharmonic Green 
function under Dirichlet boundary conditions in arbitrary C 4,7 -smooth do- 
mains. Maximum principles do not exist for fourth order elliptic equations 
and the Green function may change sign. It prevents using a Harnack in- 
equality as for second order problems and hence complicates the derivation of 
optimal estimates. The present estimate is obtained by an asymptotic analysis. 
The estimate shows that this Green function is positive near the singularity 
and that a possible negative part is small in the sense that it is bounded by 
the product of the squared distances to the boundary. 



1. Introduction 

It is well-known that the Green function G(x,y) for second order elliptic equa- 
tions on bounded domains can be estimated from above and from below by positive 
multiples of the same (positive) function, which is explicitly given in terms of the 
distances to the boundary, d(x), d(y), and the distance \x — y\. See for example [3J- 
The behaviour of the biharmonic Green function for Dirichlet boundary conditions 
should be somehow similar but will have two crucial distinctions. The singularity 
of course is of lower order, namely n — 4 instead of n — 2 with n the dimension, 
but a more serious distinction is the fact that the biharmonic Green function is 
not everywhere positive for most domains. Indeed, the few known domains with 
a biharmonic Green function of a fixed positive sign are balls in arbitrary dimen- 
sions, small perturbations of those balls and of some limacons in 2 dimensions. See 
respectively [2j [Til HTJ] and (6] [7] . The results in [7] extend and correct [13] . 

It has been observed numerically on domains with a sign changing biharmonic 
Green function that the negative part is rather small and that it is also not located 
near the singularity. The aim of this paper is to give optimal estimates from below. 
Previous results concerning smallness of the negative part have been obtained in 
[TO] for n > 3 and [4] for n = 2. With the estimates for the absolute value of that 
Green function in [5] these previously known estimates are, when n > 4, as follows: 

(1) - cd{xfd{yf < Gn{x, y) < c* \x - y\ 4 ~ n min ( 1, d ^ 2 ^f ) 

V \x-y\ j 

for all x, y G £1 C K." and where c, c* are some positive constants only depending 
on the domain. The distance of x to the boundary <9f2 is defined by 

d(x) := inf {\x — x* \ ;x* € <9£1} , 



Date: January 13, 2013. 

2000 Mathematics Subject Classification. Primary 35B51; Secondary 35J40, 35A08. 



2 



HANS-CHRISTOPH GRUNAU, FREDERIC ROBERT, AND GUIDO SWEERS 



and Gci denotes the said Green function. Let us remind the reader that the bihar- 
monic Green function Gq is such that 



u(x) 



Gn(x,y) f(y) dy 



solves 
(2) 



A 2 u = f 
= |Vw| = 



in £7, 
on d£l. 



In ([I]) the dimension is restricted to n > 4. As has been shown in [10] for n — 3, 4 
and in 4J for n — 2, the estimate from below in holds in all dimensions. An 
estimate from above has also been proved for n < 4 but the formula for that 
estimate is different from ([!]). Those estimates can be found in [5] and contain the 
function Hq in Q. 

The main result is the estimate from below in the following theorem. For the 
sake of completeness we include the estimate from above. 

Theorem 1. Let C R" (n > 2) be a bounded C An -smooth domain. Let Gq 
denote the biharmonic Green function in Q for Q). Then there exist constants 
C\ > 0, C2 > 0, depending on the domain fl, such that we have the following Green 
function estimate: 



(3) 



- 1 H n (x, y) < G a (x, y) + Cl d{x) 2 d{y) 2 < c 2 H n (x, y) 



for all x, y € £1, where 



(4) 



Hn(x,y) 



log 1 



n min <^1, 

d(x) 2 d{yf 



d{x) 2 d{yf 



\x - 

d{x) 2 - n/2 d(y) 2 - n l 2 m^ { I. 



ifn>4, 
if n = 4, 



d{x) n / 2 d{y ) n l 2 
\x - y\> 



if n= 2,3. 




Figure 1. The dashed curve gives the typical behaviour of Gq; 
the shaded area describes the band given by 

The estimate from above follows from [51 [UJ [TB] so that only the estimate from 
below has to be proved here. One should observe that G^{x,y) + c\d(x) 2 d{y) 2 > 
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was proved for some suitable c\ > in [10 . The preceding result may be con- 
sidered as an extension of the estimates in [TO] showing that close to the pole the 
positive singular behaviour of the fundamental solution can also be seen in the 
Green function. 

Due to the different behaviour of the Green function we need to distinguish 
between n > 3 and n = 2 in proving Theorem [TJ Finally we should remark that the 
lack of a maximum principle not only results in sign changing Green functions but 
also complicates the proof of these estimates in the fourth order case. The proof 
in the second order case heavily depends on the Harnack inequality which in turn 
depends on the maximum principle. 

An interesting consequence of Theorem [T] is a uniform local positivity result. 
When n > 3 this was proved in [10] . while for n — 2 we refer to [5J Theorem 6.15]. 
Here the emphasis is on the interplay between Theorem [T] and the following result. 
Moreover, we provide a proof for n = 2 which is much simpler than the previous 
one and in the same spirit as for n > 3. 

Theorem 2. Let C R" (n > 2) be a bounded C 4 n -smooth domain. Let Gq 
denote the biharmonic Green function in £1 for Q). Then there exist a constant 
rji > 0, such that 

(5) Gn(x, y) > for all x,y € f2 with \x — y\ < ra- 



2. Some auxiliary results for n > 3 



A careful inspection of the proofs in Nehari [17] and Grunau-Sweers [12] shows 
the following local estimate for the biharmonic Green function from below. 

Proposition 3. Let n > 3. Then there exists constants 5 n > and C3 > 0, which 
depend only on the dimension n, such that the following holds true. 

Assume ft C M n to be a C 4 ' 7 -smooth bounded domain and let Gq := G^2 q 
denote the Green function for the biharmonic operator under Dirichlet boundary 
conditions. If 



(6) 

then we have 



\x-y\< S n max{d(x), d(y)}, 



Ga(x,y) > < 



C3 F - V 
c 3 log ( 1 



4— n 



For the constant 5 n 



c 3 d(x)^d(y)^ 
one may achieve that S n > 0.5. 



ifn>4, 
ifn = A, 
if n = 3. 



In dimension n = 2 it seems impossible to achieve a global linear dependence of 
the radius of a ball of guaranteed positivity on the boundary distance, see Lemma 
[7]and 17 . In other words, the best result seems ^2 = ^(fi), which strictly depends 
on fi. 

Due to Proposition [3] and using the same constants S n > as there we may 
restrict ourselves in what follows to x, y such that 



(7) 



x, y G 0, x ^ y, \x - y\ > 5 n max{d(x), d(y)}. 
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Lemma 4. Suppose that n > 3 and that fl c R™ is a bounded C'" 1 -smooth domain. 
Then for each Xq G f2 t/iere exists a radius r — r Xo > and a constant C = C Xo > 
smc/i i/ia£ for all x, y G f^o.r := f2 fl B r (xo) subject to condition one has 

(8) Go(x,y)>C|x-yr"d(x) 2 %) 2 . 

Proof. We only need to discuss xo G 917 since for interior points one may choose 
r — r Xo > so small that condition |7| becomes void. 

We assume by contradiction that there exist Xk,yk G &x ,i/k — ^ ^ -^l/fc^o) 
subject to ([7]) such that 

(9) Gn(x k ,y k ) < hx k - y k \~ n d{x k ) 2 d(y k ) 2 . 

In particular we have x k — > %o, y k —> %q> d(xj~), d(y k ) — > 0, \x k — Vk \ — > 0. Without 
loss of generality we may assume that xq = and that the first unit vector e\ is 
the exterior unit normal to <9fl at xq. 




Figure 2. and subdomain Vt Xo ,i/k for xq € cT2. 



For fe large enough, we may define x k G <9fl as the closest boundary point to Xk- 
We introduce the rescaled biharmonic Green functions 

77) := - 2/fcr~ 4 G n (£fc + Nfc ~ 2/fc|£,£fc + \%k - Vk\v) 

for 

£, 77 G fife := 1 r (-ife + fl) . 

Since Xk — > 0, the exterior unit normal at <9fl converges to the first unit vector and 
so we conclude that 

f2fc — > % := {x : x\ < 0} locally uniformly for k — > 00. 

It was proved in [101 Lemma 7] that locally uniformly in H X ~H \ {(£, £); £ G 

where £* = (— £1,^2, ■ • • i£n) an d e„ is the n— dimensional volume of -Bi(O) C K™. 
We remark that this step and in particular the required uniqueness proof for Gu 
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were carried out in Grunau-Robert 10J using the assumption n > 3. The neces- 
sary modifications for n = 2 are emphasized in the proofs of Lemmata [7] and [6] 
Assumption ^ gives 

(10) G k (£ k ,T] k ) = \x k ~ y k \ n ~ 4 G n (x k ,y k ) < \\x k ~ y k \~ i d(x k ) 2 d(y k ) 2 , 

k 

where 

£fc = i Axk-Xk), Vk = -, Ayk-Xk), 

\x k -y k \ \x k -y k \ 

\Xk-Vk\ °n 

After passing to a further subsequence we hnd £, r\ € % with £ = lim^oo 
j] = lim^-^oo r] k . In view of the local smooth convergence of G k to the biharmonic 
Green function G-^ in the half space H, Boggio's formula and | — 7-7 1 = 1 there 
exists a positive constant a > such that for k large enough: 



G k (^ k ,r] k ) > ad{£h) d(rj k ) = a 



d(xfe) \ / d(y k ) 



x k -y k \J \\x k -y k \ 

4 j/, \2j/„ \2 



v\xk ~ Uk\ d(x k ) d(y k ) 



This contradicts ( 10 ) and the proof of the lemma is complete. □ 



3. Proof of the main estimate for n > 3 

Supposing that ^ holds, i.e. \x — y\ < S n max{d(cc), d(y)}, one finds that even 
(1 — S n )\x — y\ < 5 n mm{d(x), d(y)} and hence 

d(x)d(y) l~S n 
\x-y\ 2 ~ 51 

So in that case the estimate from below in Theorem [T] follows directly from Proposi- 
tion]^] Hence, we may assume from now on, again, that x, y are subject to condition 
([7]). Applying a compactness argument to 

n = [J ci xo , rxo /2 

we see that there exist positive numbers r > 0, C4 > 0, such that \x — y\ < r implies 
that Gn{x,y) > Ci\x — y\~ n d(x) 2 d(y) 2 . If \x — y\ > r we take from 5 , cf. also [9], 
that Gn(x, y) > -c 5 \x - y\~ n d(x) 2 d(y) 2 so that 

G n (x,y) + 2c 5 \x - y\~ n d{x) 2 d{yf > c 5 \x - y\~ n d{x) 2 d{y) 2 . 

Since \x — y\~ n < r~ n we end up with 

G a (x, y) + c 6 d(x) 2 d(y) 2 > c 5 \x - y\- n d(x) 2 d(y) 2 

and positive constants c§,Cq > in this case. The proof of Theorem [T] for n > 3 is 
complete. □ 
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4. Auxiliary results for n = 2 

Lemma 5. Let H = {x E R 2 : x x < 0} and let G E C 4 (H x H\ {{x, x) : x E H}) 
be a biharmonic Green function with Dirichlet boundary condition, that is 

G(x, . )A 2 ip dy = <p{x) + j (aG(x, . )B v tp ~ <pd u AG(x, . )) da 



in Jdu 

for all ip € C 4 (H) and x G Q. Moreover, we assume that G(x,y) = G{y,x) for all 
x =/= y and that a growth condition holds at infinity: 

(11) \G(x, y)\ < C (1 + \x\ 2 + \y\ 2 ) (l + (log \x\) + + (log \y\) + ) . 

Then G is uniquely determined and given by Boggio 's formula [2] • 

1 2 r\ x *-y\/\ x ~y\ v 2 — \ 
G(x,y) = G H {x,y) = — \x - y\ / dv. 

with x* — {—X\, x 2 ) ■ 

Proof. We choose some arbitrary x € H and keep it fixed in what follows. We write 

G{x,y) = G H {x,y) + H(x,y), 

where H is a regular function in H x H and A 2 H(x, . ) = in H, H(x, y) = 
d Vl H(x,y) = for y 1 = 0. According to [HI Q3], or by checking directly, with 

y* = (-2/1,2/2), 

( H(x, y) if 2/1 < 0, 

H*(x,y):=i , dH 



H(x,y*)-2 Vl ( ^- ) (x,y*)~y 2 (A y H) (x,y*) ify x >0, 



satisfies H*(x, . ) E C 4 (M 2 ) and is biharmonic on M 2 . Since G-n(x, y) satisfies ( 11 1, 
so does H{x, y). So we have 

\H{x,y)\<C x (l + \y\ 2 {\og\y\) + ). 

Using local elliptic estimates and their scaling behaviour for biharmonic functions 
satisfying Dirichlet boundary conditions on dH, that is, 

C 

\\D a H(x, ■)\\ L ac,( BRnn ) < \\ H ( X ^ ■)\\L°°(B 2R nH) 

we find that 

\V y H(x, y)\<C x (l + \y\ (log \y\) + ) , \V 2 y H(x, y)\ < C x (l + (log |y|) + ) . 

Having these estimates for H(x,y) and these first two derivatives, we have an 
estimate for H*(x, y) and may repeat the above arguments to find similar estimates 
for the derivatives of H* (x, y) and even 

|v;r(,,»)|< C , £^ . 

l + \y\ 

The maximum principle applied to the harmonic function V ' y A y H* (x, . ) shows that 
||V A mx (l + |log|fl|l) 

WVy&ytl [X,. )\\ c o {Br(0)) < O 1 + \ R \ ■ 

Letting R — > 00 yields 

VyA y H* (x,.) = 0, A y H* (x,.) = a(x) 
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with a suitable function a(. ). This shows that any WyH*(x, .) is harmonic and, 
as shown above, \7 y H*(x,y) — > as y — > oo. Hence, any X? y H*(x, . ) = and by 
Taylor's formula and observing the boundary data we conclude that 

H*(x,y) =b(x)y\ 

with a suitable function b(. ). By symmetry H(x,y) — H(y,x) and so H(x,y) — 
b(y)x\ — b(x)y\ = bx\y\, where b £ E is a suitable constant. Finally, the growth 
condition leads to b = 0, . ) = 0, and G = G-u- O 

Lemma 6 (Estimates near the boundary). Suppose that n — 2 and that C M 2 
is a bounded C^ 1 -smooth domain. Then for each xq £ dQ there exists a radius 
r — r x a > and a constant C — C Xo > such that for allx,y £ £l Xo .r '■= Q(~\B r (xo) 
one has 

(12) G n (x,y) > Cd(x)d(y) min (l, ^Mg j , 

I \x-y\ 2 ) 

Proof. We assume by contradiction that there exist Xk, yk £ Q Xa .i/k — finBi/^io) 
such that 

(13) G n (x k ,y k ) < U(x k )d(y k )min{l, d ^ d ^ 



k \ ' \x k - y k \ 2 

In particular we have x k — > xq, y k — > xq, d(x k ), d(y k ) — > 0, \x k — y k \ — > 0. Without 
loss of generality we may assume that xq = and that the first unit vector e\ is 
the exterior unit normal to dtt at xq. After possibly passing to a subsequence it is 
enough to consider one of the following two cases. 

First case: \x k - y k \ > \ Taax.{d{x k ),d(y k )} . 

This proof is as above for Lemma [4j only proving the required uniform bounds for 
the G k is slightly more involved. The arguments are sketched below in the second 
case. Thanks to Lemma [5] the convergence proof of [TU1 Lemma 7] can be extended 
to n = 2. One should observe that also the symmetry carries over to the limit. 

Second case: \x k — y k \ < | max{d(x k ), d(y k )}. 
Observe that in this case 

d(x k ) < 2d(y k ) < 4d(x k ) 

and 

\x k -y k \ < min{d(x k ),d(y k )}. 
The assumption gives that 

(14) G n (x k ,y k ) <yd(x k )d(y k ) 

k 

In this case we rescale differently, however, x k £ 9f2 denotes again the closest 
boundary point to x k . We introduce the rescaled biharmonic Green functions 

G k {i, v) : = d(x k y 2 Gn(x k + d(x k )£, x k + d(x k )n) 

for 

n k := -— (-x k + fi) . 
d(x k ) 
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Since x k — > 0, the exterior unit normal at <9£1 converges to the first unit vector and 
so we conclude that 

ilk — > H ■= {x : x\ < 0} locally uniformly for k — > oo. 

For 

6e = -77 — ^{x k ~x k ), Vk = -y — AVk-Xk) 
d(Xk) d(x k ) 



the assumption ( 14 1 transforms into 

1 2 

(15) G k (£,k,Vk) < yd k {£k)dk(Vk) < y, 

k k 

where dk ■= d( . , 9fifc). Since £k,Vk are bounded and their boundary distances are 
uniformly bounded from below by 1/2 we find after passing to a further subsequence 
that £/. — > ^ € H, rjk — > 77^ 6 H. We claim that we have local uniform convergence 
in H x H (including the diagonal) of G k to G^, since we are in dimension n = 2. 
To see this we observe first that Krasovskii-type estimates (see [131 US] and also 
Theorem 4.20]) yield at a first instance useful information only for the third 
derivatives. We have 

C 

W\ v G k {i,r])\ < 1 uniformly in k. 

Making use of 

^£dn k ,n£n k : V 6 V v G k (^r]) = 0, V 2 rj G k ^,v)=0 

and of 

VCe^ecK^ : V c V,G fc (C > 77)=0, V 2 G fe (£,77) = 
we obtain upon integration that 

|V| !7) G fe (e, ry)| < G(l + (log |£|)+ + (log |ry|) + + | log |£ - t?||) uniformly in k 
and further that 

|V^G fc (^)l < G(l + (log |^|)+ + (log|r?|) + )(l + |e| + M) uniformly in k; 
|G fc (£, 77)1 < G(l + (log + (log |r ? |)+)(l + |£| 2 + M 2 ) uniformly in k. 



Now, one may proceed further as in [101 Lemma 7]. So, ( 15 1 yields G-u{^oo,Voo) < 0, 
while Boggio's formula shows that G%(£, 77) > (even if £ = rj) since £, 77 € "H are 
interior points. □ 

Lemma 7 (Estimates in the interior). Suppose that n — 2 and i/ia£ C K 2 
is a bounded G 4 '' 1 '-smooth domain. Then for each xq G f2 t/iere exists a radius 
r = ra; > arid a constant C = G Xo > stjc/i that for allx,y £ £l Xo , r '■= H(~^B r (xo) 
one has 

(16) G n (x,y) > Cd(x)d(y) min (l, , 

Proof. Since x S fi, we have that Gq(x ,x ) > (see [T71 p. 115]): since Gq is 
continuous, there exists r, c > such that B r (xo) CC f2 and Gn(x,y) > c for all 
x,y £ B r (xo). This yields (16) since f2 is bounded. □ 
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5. Proof of the main estimate for n = 2 
Combining Lemmas [6] and [7] and applying a compactness argument to 

& = IJ tyz ,r xo /2 

we find: 

Corollary 8. Suppose that n = 2 and that ft C M 2 is a bounded C 4 '' 1 -smooth 
domain. Then there exists a radius r > and a constant C > smc/i i/iai /or aZi 
.To 6 O and /or a// x, y € Q Xo ,r := H -B r (a;o) one /ias 

(17) G n (s,g) > Cd(s)d(y) min (l, f x)d ^ ) . 

This step provides in particular a different and simpler proof for the local pos- 
itivity statement from [51 Theorem 6.15] which was proved first by Dall'Acqua, 
Meister, and Sweers |J. 

Proof of Theorem^ If \x — y\ > r we take from [5], see also [9], that Ga(x,y) > 
-c s \x - y\~ 2 d(x) 2 d(y) 2 so that 

G n (x,y) + 2c 8 \x - y\- 2 d(x) 2 d(y) 2 > c 8 \x - y\- 2 d{x) 2 d{y) 2 

> C8 d(.)d( y) min{l,f^j. 

Since |ac — y| ~ 2 < r~ 2 we end up with 

G n (x, y) + c 9 d(x) 2 d(y) 2 > e 8 d(x)d(y) min (l, d ^ d ^ \ . 

and positive constants c 8 ,cg > in this case. The proof of Theorem [T] is now 
complete also for n = 2 using Corollary [8] □ 

6. Proof of Theorem [2] 

This theorem was proved in [10] when n > 3 and in [9, Theorem 6.15] when n = 2. 
The latter proof is quite involved and based on the extensive use of conformal maps 
and explicit Green functions in certain limagons. We provide here an alternate proof 
which uses the same techniques for n — 2 as for n > 3. 

Case I: d(x)d(y) < \x — y\ 2 . For this situation we have 



(18) H n (x,y) = 



x~y\- n d{xYd(yY iin^A, 

d(x) 2 d(y) 2 \ 
log [ 1 + ; V f; I if n = 4. 

\x-y\ 4 



Then there is c > such that we find c 2 1 Hq(x, y) > c\d{x) 2 d(y) 2 for \x — y\ < c. 



Case II: d(x)d(y) > \x — y\ 2 . Now we have 



(19) H n (x,y) 



|a; — 2/1 ifn>4, 
, / d(x) 2 d(y) 2 \ 

4 d(x) 2 - n/2 d{y) 2 - n/2 ifn = 2,3. 
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Since d(x)d(y) is bounded on fi one finds for n > 4 the existence of c > such that 
Ha(x,y) > c\d{x) 2 d(y) 2 for \x — y\ < c. For dimension n — 3 the argument is 
more subtle. We fix 6 n as in Proposition M Taking d{x)d(y) < e 2 ^ for e > 
but sufficiently small it follows that c^ 1 Ha(x,y) > c\d{x) 2 d{y) 2 . It remains to 
consider d(x)d{y) > ^^s 2 . Assume first that \x — y\ < ^fe and d(x) < |e. Then 
d(x)d(y) < 1+ 2 n £ 2 and we are in the situation just considered. So we are left with 
\x — y\ < 4f-£ and > \e. Then we may apply Theorem 1 of [12] . see also 
Proposition [3J to find that for |a; — y| < |# n £ < <5„ max(d(x), it follows that 
Gq(x, y) > 0. For dimension n = 2 the result follows directly from Corollary [8j 
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